NEW  YORK  UNIVERSITY 

INaTlTUTE  OF  MATHEr,!/\TICAL  SCIENCES 

AFCRC-TN-56-798  LIBRARY 

ASTIA  DOCUMENT  No.  AD  11 01 39  15  ^»k^  Pbce,  New  York  3,  H% 

«.£  ET  PR/I^ 
^f.T\r%     NEW   YORK   UNIVERSITY 

^   I  In   iJL    ^     Institute  of  Mathematical  Sciences  ^ 

"^      v--/        <      Division  of  Electromagnetic  Research  3 


RESEARCH   REPORT   No.    EM-99 


Excitation  of  Surface  Waves 


B.  FRIEDMAN  and  W.  ELWYN  WILLIAMS 


CONTRACT  No.  AF  19(604)1717 
OCTOBER    1956 


J 


NEW  YORK  UMIVERSITT 

Institute  o.f  Mathematical  Sciences 
Division  of  Electromagnetic  Research 

Research  Report  No,   01-99 


EXCITATION  OF  SIRFAGE  WAVES 

i  ; 

B.  Priedmaii  and  W,  Elwyn  Williams 


B.  Friedman 


W.  Elwyn  Williams 


Morris  Kline 
Project  Director 
October,  1956 


The  research  reported  in  this  dociment  has  been  made  possible 
through  the  support  and  sponsorship  extended  by  the  Air  Force 
Cambridge  Research  Center,  under  Contract  No.  AF19(60U)1717. 
It  is  published  for  technical  information  only  and  does  not 
necessarily  represent  recommendations  or  conclusions  of  the 
sponsoring  agency. 

New  York,  1956 


Table  of  Contents 


Page 


lo   Introduction  1 

2.  Outline  of  results  3 

3,  Mathematical  formulation  8 
Ut       Interpretation  of  solution  11 

5.  Asymptotic  behavior  of  the  field  13 

6.  Comparison  of  power  16 

References  26 

Distribution  List 


Abstract 


It  is  shown  how  to  locate  a  dipole  source  above  a  dielectric 
surface  so  as  to  either  (a)  produce  as  pure  a  surface  wave  as  possible 
or  (b)  majdmize  the  amount  of  energy  carried  by  the  surface  wave. 
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1,       Introduction 

Since  the  appearance  of  iSommerf eld's  classical  paper L  J  on  the  field  produced 
by  a  dipole  above  a  plane  dielectric  earth,  a  considerable  amount  of  attention  has 
been  devoted  to  the  problem  of  the  excitation  of  surface  waves  over  dielectric  sur- 
faces.    The  two  questions  which  are  of  practical  importance  are  these:     l)  to  deter- 
mine under  what  conditions  a  surface  wave  exists  and  2),  if  it  exists,  to  determine 
the  type  of  source  which  is  the  most  efficient  generator  of  the  surface  wave. 

Though  it  is  now  generally  accepted  that  the  surface  wave  due  to  a  dipole  abo-ve 
a  dielectric  earth  is  nowhere  predominant  on  the  earth's  surface,  the  question  as 
to  whether  there  exists  a  source  distribution  which  produces  a  surface  wave  pre- 
dominant in  some  region  on  the  earth's  surface  still  remains  unanswered.     To  answer 
this,   an  investigation  of  the  propagation  above  a  plane  earth  when  the  source  dis- 
tribution is  arbitrary  is  necessary.     This  problem  has  been  considered  both  by 
Kay"-  ■'   and  one  of  the  authors'-  J   and  it  has  been  shown  that  the  surface  wave  and 
the  radiation  field  t3rms  may  be   obtained  from  a  knowledge  of  the  far-field  radia- 
tion pattern  of  the  original  source  distribution.     From  these  results  it  is  theoretic- 
ally possible  to  combine   sources  with  different  far-fielc  patterns  in  order  to  deter- 
mine the  most  suitable  source  distribution  for  generating  surface  waves,     Since  this 
approach  requires  a  large  amount  of  tedious  nujnerical  work,  it  seems  desirable  to 
attack  the  problem  in  a  more  analytic  manner. 

This  haj  been  done  by  Cullen  in  a  recent  paper'--'.     There  he  considers  the  ex- 
citation of  surface  waves  by  a  horizontal  slot  placed  above  a  plane  surface  on  which 
an  impedance  type  of  boundary  condition  holds.     For  a  reactive  type  boundary  condition 
it  is  known  that  a  surface  wave  exists.     Cullen  discusses  two  different  criteria  for 
determining  that  height  of  source  above  the  ground  plane  which  is  most  efficient  in 
exciting  the  surface  wave.     By  the  first  criterion  the  most  efficient  soxnrce 
is  that  which  produces  an  almost  pure  surface  wave   (this  concept  will  be  explain- 
ed later)  near  the  guiding  sxirface.  By  the  second  criterion  the  most  efficient   source 
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is  the  one  which  produces  a  field  such  that  the  ratio  of  energy  per  unit  area  of 
wavefront  in  the  surface  wave  to  the  sum  of  energy  per  unit  area  of  wavefront  in 
the  radiated  field  and  the  siirface  wave  is  a  maxiinuin.     Cullen  obtains  values  of 
the  optimum  source  height  for  each  of  these  criteria.     However,  when  discussing 
the  second  criterion,   he  finds  it  necessary  to  use  nxunerlcal  methods  aril  thus  ob- 
tains the  optimum  source  height  only  for  one  particular  value  of  the  surface  im- 
pedance. 

The  present  work  was  suggested  by  Cullen 's  paper  but  it  is  more  general  in  that 

both  two- and  three-dimenaional  problems  are  investigated  and  a  more  general  type  of 
guiding  surface  is  employed.     The   surface  is  assumed  to  be  a  slab  of  dielectric 
resting  on  a  ground  plane  on  which  an  impedance -type  boundary  condition  holds. 
The  dielectric  need  not  be  homogeneou8|  in  fact,   it  may  consist  of  layers  of  different 
materials^  but  we  assume  that  the  permittivity  of  the  dielectric  varies  only  in  the 
direction  perpendicular  to  the  ground  plane.     When  investigating  the  field  produced 
by  a  source  above  this  guiding  surface,  we  find  it  desirable  to  expand  the  source 
into  plane  waves  of  all     wavelengths  in  a  direction  perpendicular  to  the  ground 
plane  and  not  into  plane  waves  in  a  direction  parallel  to  the  groxind  plane,  as  Cullen 
does.     With  our  type  of  expansion  we  find  that  the  only  property  of  the  guiding  sur- 
face that  we  need  is  its  impedance  Z(k)  to  an  incoming  normally  incident  plane  wave 
of  wave  number  k.     The  optimum  height  h^  above  the  surface  at  which  a  slot  must  be 
placed  to  provide  a  nearly  pure  surface  wave  is  then  determined  to  be 

L 

o 

where  k     is  the  wave  number,   Z     the  characteristic  impedance  of  free  space, and  where 
0  '     o  "^  ' 

we  put 

z(k)     »     K(kj  ♦  i  X(k), 

80  that  A(k)  is  the  reactance  of  the  guiding  surface.     This  result  generalizes  that 
obtained  by  Cullen. 


-  3  - 

We  consider  also  a  second  criterion  for  launching  efficiency  which  differs  from 
that  proposed  by  Cullen.     The  criterion  is  as  follows:     That  height  h^  of  source  is 
most  efficient  for  which  the  energy  in  the  radiated  field  per  unit  area  of  wavefront 
is  a  minimum.     If  Z(k)  is  purely  reactive  we  find  that  h^  is  given  by  the  following 
formula:  / 

In  the  case  considered  by  Cullen  these  two  criteria  give  practically  the  same  re- 
sults for  h„  as  those  he  obtained, 

Ihe  formula  for  finding  the  most  efficient  launching  heigt^t  h^  may  be  given 
the  following  physical  interpretation:     Consider  a  plane  wave  with  wave  number  K 
travelling  toward  the  interface  in  a  direction  making  a  U5     angle  with  the  vertical} 
then  h„  is  that  height  at  which  the  total  field  produced  by  this  wave  and  the  re- 
flected wave  is  a  minimum. 


2.       Outline  of  results 

Suppose  that  the  dielectric  material  extends  from  z  ■  -z     to  z  «  0  and  that 
there  is  free  space  for  z  >  0.     We  suppose  that  the  exciting  source  is  a  horizontal 
magnetic  slot  parallel  to  the  y-axis  and  passing  through  the  point  z  ■  h,  x  »  a. 
The  field  quantities  can  be  expressed  in  terms  of  H     as  follows: 

H     -     (0,Hy,0) 

E     .     (icoer-^  (  ^  ,0,  -  1^  )  Hy  . 

There  are  two  methods  for  obtaining  H  .  One,  which  Cullen  uses,  is  to  expand  the 
source  in  terms  of  waves  travelling  along  the  x-axis.  The  other,  which  is  used 
in  this  paper,  is  to  expand  the  source  in  terms  of  waves  which  come  in  and  are 
reflected  along  the  z-axis. 

Suppose  that  a  wave  e     comes  in  along  the  z-axis.  Because  of  the  presence 

of  the  dielectric,  this  wave  will  produce  a  reflected  wave  p(k)  eikz  for  z  >  0, 
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Here  k  will  be  any  reaj.  nximber  from  0  to   oo.     However,  if  the   source  is  oscillating 
with  frequency  co.thus  producing  waves  whose  wave  nuiiiber  in  free  space  is  k  ,   then 
the   incoming  wave  will  be 

e-^^^exp[iy^|x.al] 

and  the  reflected  wave  will  be 

p(k)  e^^*  expfi  Vk^  |x.a|  J    . 
Notice  that  this  may  be  interpreted   as  a  wave   coming  in  at  an  angle  9  to  the  x-axis 
and  then  reflected  at  the  x-axic   (see  diagram   below).      Here  sin  ©  '   ^/^q' 


When  W8  take  into  account  the  fact  that  the  source  producas  plane  waves 
travelling  in  all  directions  we  find   (see   (12')     that  H    will  contain  a  tera 


proportional  to 

,00 


r.    /JT7  ,        ,1  ik(z-h) 
expl  i  /k^-k     |x-a|Je 


dk 


..-^—T 


^- 


00  r    /T~7  ,        ,1  ik(z-*-h) 

/^  expi/k-k      |x-a|   e 

-      p(^)  L_^ =! dk  . 

'Loo 


The  first  integral  is  just  the  field  that  would  be  produced  by  the  source  in  free 
space  if  the  dielectric  were  not  present^  and  the  second  integral  is  the  reflected 
field  that  is  produced  by  the  dielectric.  The  sum  of  these  two  integrals  will 
give  the  radiation  field* 
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Let  UB  consider  p(k).     There  may  be   a  complex  value  of  k  with  positive 
imaginary  part,  call  it  k  ,  such  that  p(k.)  is  infinite.     This  value  of  k,  will 
produce  a  surface  wave  and  we  shall  have  in  H     (see  (i2))  a  term  proportional  to 

exp     i /k^-k^  |x-a|     e     "^    e  ,         z  >  0. 


Notice  that,  because  the  imaginary  part  of  k,   is  positive,  this  expression  represents 
a  wave  which  travels  along  the  x-axis  at  a  speed  faster  than  light  and  is  exponentially 
attenuated  in  the  z-direction.     This  behavior  is  typical  of  surface  waves. 

The  value  of  the  reflection  coefficient  can  be  found  from  the  impedance  of  the 

— ikx 

dielectric  to  an  incoming  wave  e         .We  have,  from  (17), 

Z(k  sin©; 
sin  ©  + 


Z^ 
p(k  sine)  -  ° 


0  z(k  sin©) 

sin  ©  - 


^0 


From  (iii)  the  radiation  field  at  large  distances  R  from  the  source  and  in  the 
direction  ©  to  the  x-axis  is  proportional  to 

^^0^   r--ik  ksin©  ik  ksin©- 

e 


O  1 IK    KSXnW  XK.     "l.S-LUW-1 

e  +  p(k  sin©)e 

Tr      L  °  J 


Using  tlds  we  find  (see   (25))  that  the  radiated  power  per  unit  area  of  wavefront 

is  proportional  to 
n 

/.  2ik  ksin©     2 

I  1  I-  p(k^sin©;  e      °  d©. 

o 
This  formula  may  be  given  a  simple  physical  interpretation.     The  field  produced 

by  the   source  may  be  decomposed  into  plane  waves  travelling  in  all  directions.     The 

plane  waves  which  travel  in  the  direction  -©  will  be  reflected  back  to  the  height  h 

2ik  ksin© 
of  the  source  with  a  reflection  factor  p(k  sin©;e  .The  total  field  at  a 

distance  R  will  then  be  the  sum  of  the  direct  plane  wave  whose  amplitude  is  A,   say, 

2ik  ksin© 
and  the  reflected  plane  wave  with  amplitude  Ap(k  sin6)e  .     The  total  energy 

power  in  the  wavefront  will  be  the  integral  of  the  absolute   square  of  the  value  of 

the  field,  and  in  this  way  we  get    the  equation   (-m-)© 
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Put  p  -    |p|  e        '^where    |o|   and  y  are  real.     Using  this,  we  find  that  the 
radiation  pattern  F(&)  is  given  by 

,   -ik  hsine  i(K  hsinS+Zy) 

F(e)  -     e       °  +    Ipl  e       ° 


1  +    |p|     +  2|p(   cos   2(lCQh  sin9  *  yU 


-.1/2 


Note  that  the  value  of  the  reflection  coefficient  can  be  obtained  from  the 

impedance  of  the  dielectric  to  an  incoming  wave  e  ,     We  have,  from  (17), 

;i(K  sin9) 
sin  ©  ♦ 


^o 


P(^'^"®^  -  Z(K  sin6) 


sin  e 2 

In  case  Z  is  a  pure  reactance,  this  formula  shows  that   |p|   ■  1,  and  a  knowledge  of 
the  far- field  pattern  F(©)  enables  us  to  obtain  y  *nd  thus  also  p(K  sin© J.     Theoreti- 
cally, this  information  about  p  permits  us  to  obtain  the  values  of  6  for  which 
p(k  sin©)  ■   00,  that  is,  for  which  surface  waves  exist.     This  agrees  with  the  re- 
sults in    [2]   and   [3],     However,  in  practice  it  is  very  difficult  to  find  the  surface 
wave  this  wsty  because  the  far-field  pattern  gives  the  values  of  p(k  sin©)  for  real 
angles  6  whereas  a  surface  wave  corresponds  to  complex  angles  ©. 

A  glance  at  the  Fibres  will  show  that,   at  least  in  the  case  where  Z  is  a  con- 
stant reactance,  the  radiated  power  per  unit  area  of  waveffont  is  approximately  pro- 

2 
portional  to  F(©)     evaluated  at  ©  «  n/U.     Using  this  approxifiHtion,  we  may  deteraine 

the  height  h„  for  which  the  radiated  power  per  unit  area  of  wavefront  is  a  miniMiB, 

We  find 

k^  hj  sin  J  ♦  Y     -     5  , 

or 

o 
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-i            X 
bince  ir  the  case  of  a  pure  reactance^  y  "  tan         ,  we  get    (see   (28)) 

L^  sin  p 


^.f  cct-^(yr   ^). 


Going  back  to  the  physical  interpretation  of  (»),  we  can  interpret  h^  as 
follows:  The  height  of  source  hp  at  which  the  radiated  power  per  unit  area  of 
wavefront  is  a  mLnimum  is  approximately  the  height  at  which  there  is  complete 
interference,  i.e.,  zero  amplitude  between  the  direct  and  reflected  waves  travelling 
in  the  U5  -direct J on  to  the  interface* 

For  experimental  purposes  it  is  useful  to  be  able  to  produce  a  field  which  is 
as  pure  a  surface  wave  as  possible,  bince  the  surface  wave  attenuates  exponentially 
in  a  direction  perpendicular  to  the  interface,  we  consider  the  f  ield  at  a  point  on 
the  interface  at  a  distance  x  from  the  source.  From  (l2)  the  H  component  produced 
by  the  surface  wave  is 

n — 

i  VK  -X.  X 

1    e    °  ^ 

or         u(o,\,;u(h,x.  ), 

\/k  -  X, 
*   o        1 

whereas  from  (l8j  the  H  component  produced  by  the  radiation  field  is 

V  n    mr   — - — yg        u^(o)u  (h), 

-5/2 
plus  terms  of  order    |x|     '    »     By  proper  choice  of  the  launching  height,  the  above 

-3/2 

term  in  the  radiation  field  of  order  |x|  '  can  be  made  zero, and  therefore  the 

total  field,  which  is  the  sum  of  the  surface  wave  and  the  radiation  field,  becomes, 

at  least  at  large  distances  from  the  source,  very  nearly  a  pure  surface  wave.  The 

height  at  which  this  term  in  the  radiation  field  vanishes  is  given  by  (19)  as 

iZ 
^  "  k  Z(0)  • 

0 
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3»       Mathematical  formulation 

Consider  a  ground  plane  at  z  ■  -z     and  suppose  that  the  permittivity  of  the 

medium  above  it  depends  on  the  z-coordinate  only.  For  example,  the  medium  might 

be  a  dielectric  slab  from  z  =  -z     to  z  =  0  and  then  free  space  from  z  »■  0  to  z  =oo: 

o 

or  the  medium  might  be  composed  of  layers  of  dielectric  of  arbitrary  thickness  and 
permittivity  from  z  =-z     to  z  »  0  and  free  space  above;  or  the  permittivity  of  the 
medium  might   vary  continuously  with  height  from  z  =  -z     to   z  =  0  and  again  have  free 
space  above.    In  all  cases  we  shall  assiane  that  "the  medium  between  z  =  0  and  z  =oo  is 
free  space  and  that  it  contains  a  source  of  electromagnetic  radifition.     Finally,  we 
assume  that  the  electromagnetic  field  satisfies  an  impedance-type  boundary  condition 
at  the  ground  plane  z  «  -z   . 

We  shall  treat  both  a  two-  and  a  three-dimensional  problem  simultaneously  with 


the  help  of  the  following  coordinate  S3rstems: 


^Z 


In  the  first  diagram  above,  the  soirrce  of  the  electromagnetic  field  is  a  magnetic 
line  dipole  placed  at  z  "  h,   x  ■»  a  and  extended  parallel  to  the  y-axis.   Because  of  the 
geometry,  the  field  produced  will  be  independent  of  y.  We  shall  call  this  the  two- 
dimensional  problem  and  denote  it  by  Case  I. 

In  the  second  diagram  above,   the  source  is  a  circular  magnetic  dipole  located  at 
z  =  h,   r  =  a.  Because  of  the  geometry,  the  field  produced  will  be  independent  of  0. 
We  shall  call  this  the  three-dimensional  problem  and  denote  it  by  Case  II.     In  both 
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The  simple-harmonic  form  for  Maxwell's  equations  is 

curl  H     -     -i  CO  e  E 
curl  E     ■     i  »  |i  H 

where  E,  H  are  the  electric  and  magnetic  field  intensities  and  e  and  ^,  are  the  per- 
mittivity and  magnetic  permeability.  Note  that  by  assumption  6  is  a  function 
of  z  and  n.  is  a  constant. 

In  Case  I,  because  the  source  is  a  magnetic  line  dipole,  the  field  is  completely 
determined  by  the  value  of  H  «  We  have 

H   -  E   -  H   -  0 
X     y      z 

and 

1  ^"v  1  ^" 


The  field  component  H     satisfies  the  equation 

1/2      '  / 

Put  H     »  6  '     H     and  (1)  becomes 

y  y 

2    '  2    • 

3  H         9  H  «    ,  •» 

(2)  — ^  +  —J.  ♦  co'e  n  H     »     i-    6(z-h)6(x-a) 
3x^         9z  y        y?      ■ 

where 

/,N  2    '  2  1  d^e       3   /d6,2    I 

(3)  CO  e  |i    -    w  eti  +  ^^    — 7  '  7  ^^'     (* 

In  Case  Il^because  of  our  assumption  about  the  source^ the  field  will  be  complete- 
ly determined  by  the  value  of  IL»     We  have 

H,     -     H^     -     Eg     -     0 
and 


-  iO 


(-ia)e)"V^  If  ^"""e^' 


The   field  component  H     will  satisfy  the  equation 


ih) 


2  2 


38         ar 


5(z-h)5(r-a) 
— 7^ 


,V2, 


Put     tL  -  c       H_  and  the  following  equation  is  obtained: 

2    '  2    ' 

^  "e       i  ''"&        /2    '  l\, 


(3) 


■Hg        d^Hg        .    9H, 


5(z-h)5(r-a) 

172— 


2nr  e 


where  e     is  defined  by  (3). 
We  set 


2    »  2  ,    , 


k^  -  q(z) 


where  e     is  the  permittivity  in  free  space.     Note  that  by  our  assiiraption  q(z)  -  0  if 

°  ^2 

z  >  0,     If  we  denote  the  operator  q(z)  -  — «•    by  L,   the   formal  solution  of  equation 


(2)  is 
(6) 


a"? 


•    -1/2 
H  -  e  ' 

y 


exp 


(iy^  (X- 


2iV^ 


=L  6(z-h) 


and  the  corresponding  formal  solution  of  {5)   is 


(7) 


.-■»-/2 


^e'  "uT 


''i(>^V^*)  "i""  V^^^O  ^^'"^^' 


r  >  a, 


These  results  will  be  interpreted  in  the  following  section. 
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U.       Interpretation  of  solution 

To  interpret  (6)  and   (7)  is  is  necessary  to  express  6(8-h)  in  terms  of  the 
eigenl'unctions  of  the  operator  L. 

Consider  the  eigenvalue  problem 

(8)  Lu(z)  -     -  — H-    *  q(z)u     ■  KM 

dz 

where  u(z)   satisfies  the  given  impedance-type  boundary  condition  at  z  »  -Z-»     If 

q(z)  and  the  impedance-type  boundary  condition  at  z  =  -z     are  real,  there  may  exist 

negative  values  of  K  such  that   (8)  has  a  non-zero  solution  satisfying  the  given 

boundary  condition  at  z  ■  ~^a*     ^^  ^(^^   (o^  'the  impedance-type  boundary  condition) 

is  not  real,  there  may  exist  complex  values  of  A.  with  negative  real  parts   such  that 

(6)  has  a  non-zero  solution  satisfying  the  given  boundary  conditionat  z  «  -z   .     In 

either  case  we  shall  denote  these  values  of  \  by  K.  and  the  corresponding  normalized 

solution  of  (8)  by  u(z,X..)«     Note  that  u(z,\,)  vanishes  exponentially  as  z  -s>oo. 

J  J 

Every  positive  value  of  K  will  be  in  the  continuous  spectrum  of  the  operator  L» 
Put  X  "  k  and  denote  by  u(z,k)  that  solution  of  (8)  which  satisfies  the  boundary 
condition  at  z  »  -z  ana  vrtiich  as  z  ->oo  has  the  asymptotic  form 

(9)  u(z,k)  -  e-^^^   +  p(K)e^*'^  . 

— ikz 
Here  p(k)  is  the  reflection  coefficient  for  the  incoming  plane  wave  e     incident 

on  the  dielectric  at  z  -  0.  Since  q(z)  =0  for  z  >  0,  formula  (9)  is  an  equality 

for  z  >  0«  We  haveL^J 

(10)  5(z-h)  -  73  u(z,A..)u(h,A..)  +  ^  j  u(z,k)  uThTFJ"  die  . 

J  -b 

Here  u(z,k)  is  that  solution  of  (8)  satisfying  the  boundary  condition  at  z  =  -z 
which  for  z  >  0  has  the  form 
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(il)  u(z,kj     -     e         +  p(-ic)e  . 

If  q(z)  and  the  impedance-type  boundary  condition  at  z  ■  -z     are  real,  the  value 
of  p(-k)  is  the  complex  conjugate  of  the  value  of  p(k)  and  hence  u(z,k)  is  the 
complex  conjugate  of  u(z,k).     However,   in  all  cases,    comparison  of    (9)   and   (ll) 
shows  that  p(-k)  -  p(k)"   . 

If  the  expression  (iO)  is  substituted  in  equations   (6)  and  (V)   and  il  we  use 
the  fact  that  u  satisfies  (8),  we  get  the  following: 

,        ^-1/2  expli/k^   tx-a|] 

"y  ■  -H-    f  -    ,,  -  u(z,^^;u(h,^j) 

'   o         j 


(12) 


,-V2 

♦  -1 — ■ — 
Uni 


I       /2     2  I 

exp  i  /K^-k      |x-a| 

/^ — :? 

yk     -  k 


u(z,kj  u(h,k;  dkj 


«C 


.       e-^/2 


^    ^    J,(^/^a^Hii)^/?Xr)u(z,..Mh,^.) 


(13) 


,-1/2       '^ 


*  "BH" 


j   Jj^  ( /V?  a  ]    H^-*-^  (/k^-k^  r  j  u(z,k)  u{h,k)  dk  . 


The  first  term  on  the  r ight-hand  side  in  equation  (l2)  or  equation   (l3)  ie 
the  surface  wave,  whereas  the  integral  term  in  these  equations  is  the  radiation 
field.     These   formulas  show  that  the  existence  of  a  surface  wave  depends  on  the 
existence  of  eigenvalues  of  (8)  with  negative  real  parte.     Thus  existence  depends 
only  on  the  value  of  q(z)  and  the  impedance-type  boundary  condition  at  the  ground 
plane,  that  is,  only  on  the  electromagnetic  properties  of  the  medium  in  the   z- 
direction  and  not  on  the  number  of  dimensions. 
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5»       Aeyinptotlc  behavior  of  the  field 

We   shall  obtain  the  asymptotic  behavior  of  the  integrals  in  (i2)  and  (13) 
and  from  these  obtain  the  asymptotic  value  of  the  electromagnetic  fieic  conponentSi 

Using  (9)  and  (ll),  we  see  that 


Putting  this  expression  into  (l2)  and  (l3),  we  find  that  the  integral  term  becomes 


(12') 


.00 


/_oo 


exp  i  Vk^-K      |x-a| 


7 — rr 


♦  ik(z-h)  -    ,'  ik(z+h) 

e  +  p(Kje     ^       ' 


dk 


in  Case  I  and 


(13') 


pi(^^«)  4'(V^  -)  E*''^'''^  *  P(K)ei'^(-^)] 


dk 


in  Case  II. 

Ihe  method  of  stationary  phase  will  be  used  to  obtain  asymptotic  expressions 
for  these  integrals.     In  Case  I  we  put 


X  -  a  «  R  cos  ©, 


z  "  R  sin  0 


but  in  Case  II  we  put 


r  -  R  cos  ©, 


z  -  R  sin  9. 


For  large  values  of  rt,  the  integral  in  (l2J  is  asymptotically 


(lU) 


.1/2       ^(''o"-  l^    ~ 
6  e 


V^       y^^ 


■ik  h  sine  *^o^  ^^^'^ 

e  +  p(k  sin©)e 


] 


and  the  integral  in  (13^  is  asymptotically 
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■1/2    ^("o^-r' 

C  6 


o 


-ik  h  sine  +ik  h  sin© 

e  +  p(k  sin0)e       ° 


The  expressions   (lii)  and   (i5)  are  essentially  the  radiation  field  parts  of 
H     and  H^  respectively.  We  now  consider  the  value  of  the  radiation  field  on  the 
surface  of  the  slab,  at  large  distances  from  the  source,  that  is,   at  z  =  0  but 
R  »  0.   In  this  case  €  approaches  zero  and  we  can  show  that  the  expressicm  in 
brackets  in   (lii)  and   {!$)   approaches  zero.  This  indicates  that  at  z  =  0  the  asymp- 
totic expression  for  the  radiation  field  is  no  longer  given  by  (lU)   and   (l5). 

We  shall  obtain  the  correct  asymptotic  fonriula  for  the  radiation  field  on  the 
surface  of  the  slab  for  R  large  by  evaluating  the  integrals  in   (12)  and  (13 )  in 
a  differait  manner.   However,  first  we  evaluate  p(k  sin©). 

At  the  interface  z  ■  0,  between  the  dielectric  slab  and  free  space,  t'r.s  ratio 

E  /H     or  F  /H^  Tnust  be  continuous.     Denote  this  ratio  by  Z,(k)jthen 
X     y  r     ©  -  » 

3H 


or 


(^^^  u  dz 


(^^^     ar  /  «y    ■    ^(*^^» 


1  du  ,    .„/...  o 


-  icoez(k)  -  —^    Z(k), 

2»0  O 


where  Z  is  the  characteristic  impedance  of  free  space.  Note  that  Z(k)  is  the 
c 

-ikz 
impedance  at  the  interface  to  an  incoming  wave  e    *  For  convenience,  we  shall 

sometimes  write  this  impedance  as  Z.  For  z  >  0  we  have 

,  ,  >      -ikz    ,,  V  ikz 
u(z,k)  -  e     +  p(k)  e    ; 

then  by  (lb) 
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or 


1-p 
i*p 

kZ^     > 

k         Z 

p    » 

0          o 

k       z 

O            0 

r  iim-Lxjf ) 

p(k  sing)     . 
o 

Z(k  sine) 

(17) 

0 

Z(k  Sin©) 

cf-i  J-.    Q    _            ^ 

0 

From  this  formula  it  is  clear  that  if  Z(0)  f  0,  the  value  of  p(k  sin©) 

approaches  -1  as  ©  approaches  zero  and  therefore  the  bracketed  expression  in 

(iU)  and  (15)  approaches  zero.  Also,  from  (l7), 

2Z 
p(k^sin9)  .  -1  -  ^^  sin  © 

for  sma].l  values  of  ©,     Then 

-ik^z  sin©  ik  z  sin© 

u(z,k^sin6)     «>     e  +  p{K  sin©)e     ° 

0  0 

Z. 


^sin© 

k 
o 


\jm  *  'v] 


We  shall  denote  the  bracketed  expression  in  this  formula  by  u-,(z). 

Let  us  return  to  the  integrals  in  (i2)  and  (l3)  and  consider  their  asymptotic 
behavior  for  large  values  of  x  or  r.  When  we  make  the  change  of  variables 


2     2     2 

k   -  k   -  B 
0    ^ 


we  haYS 

dk         d^ 
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therefore  the  resulting  integrand  becomes  infinite  at  p  -  K^,  that  is,  at  k  =  0. 
This  indicates  that  the  asymptotic  behavior  of  the  integral?  in  (i2)  and  (i3)  can 
be  obtained  by  expanding  the  integrand  in  powers  of  K,  In  this  way  we  find  that 
the  integral  in  Case  I  is 


(18) 


/F  W-  \   °   ,3/^  u,(z)u^(h; 


k^|x-a 


and  in  Case  II  it  is 

o 
If  we  choose  u-(h)  such  that  u,(h)  «  0,  then  clearly  both  expressions  for  the 
radiation  field  will  vanish.  This  means  that  to  obtain  as  pure  a  surface  wave  as 
possible  on  the  interface,  we  must  require  the  height  h  of  the  source  to  be  such 
that 

(19)         k  h  -  — ^  . 

Z(0) 

Sliice  the  source  is  always  above  the  interface,  that  is,  h  >  0,  formula  (19) 
shows  that  a  surface  wave  will  exist  only  if  Z  is  a  positive  reactance.  In  case 

Z  is  complex,  formula  (l9)  gives  a  complex  value  for  h,  which  means  that  at  no 

real  height  will  u, (h)  -  0,  However,  if  we  put  Z  ■  R  +  iX,  a  surface  wave  can 

exist  only  if  X  >  0  and  we  see  that  |\i,(h)|  will  be  a  minimum  at  that  height  h 

for  which  .. 

(20)  \^  '  *  my  ' 


6»   Comparison  of  power 

In  this  section  we  compare  the  radiated  power  to  the  power  in  the  surface  wave. 
For  simplicity  we  shall  restrict  ourselves  to  Case  I,  Similar  results  can  be  obtained 
for  Case  II. 
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From  (12)  the  amplitude  ol"  one  siirface  wave  term  in  H     is 

iyi^|x-a| 

'  o         1 

Note  that 


(21)  (2i)~^  ^ u(z,X^)  n{h,\^) 


1  /P     ^  AT  ^ 
(22)  u(z,>^)     -     {2±s/K[r'''  e         ^       ,  z>0. 


The  x-component  of  the  average  power  ie 


1  -1   r  ^«y      * 

(23)  |Re(i(oe)  ^  ^f    ^  '^^ 


where  H     is  the  complex  conjugate  of  H  • 

5ince  power  orthogonality  between  surface  wave  and  radiated  field  occurs  only 
in  case  >_.   is  real,   we  shall  assume  henceforth  that  X_   is  real.     Then,   applying 
(23)  to  (2l),  we  find  that  the  average  power  in  the  surface  wave  is 


(2U)  _-L_    e-2Hh 

if  K^     -     -p^. 

To  find  tha  radiated  power  Jar  unit  area  of  wevefront  we  must  evaluate 

n/2 
^    KB  H*  E^  dO 


or,  since 

-1  ^"y 


we  must  evaluate 


-  IS  - 


7    Redcoe)-^    j         ^    h; 


d9. 


Using  (lU)  multiplied  by  e  ^     for  the   amplitude  of  the  radiated  field  in  H  , 
find  that  the  radiated  power  per  unit  width  of  wavefrcnt  is 


we 


l6nK 


-1      (^^"^  I ik  h  sinQ 


+  p(k  sin©)e 


ik  h  sin9 
o 


] 


(25) 


r 


ik   h  sin©  "^^r^  3in9- 

j     °  +  p(-k^sin9)e  I    d9. 


'] 


If  we  assume  Z  is  purely  reactive  and  constant  for  all  K,  we   shall  get  Cullen's 
results.     Put  Z  -  iaZ  ;  then  from  (17) 


p(k  sin9) 


i  8in9-a 
i   sin9+a 


and   (25)  becomes 


(l6n<oeK   )"■'■  I(h,h) 


where 


Kh 


.''/^ ik  h  sin9  .      .    -     ik  h  sin9 

a-i  sin9         o 
e 


a+i  sin9 


] 


(26) 


[ 


ik  h  sin9     ^^^  g.^e     -iK^ 


-ik  h  sin9-| 


a-i  sin9 


d9. 


It  is  easy  to  show  that 


ly 


(^^'^  r        (  1  1  A    -^         -^^  sine(h+h')-i 


dO 


-  «   •^oj^o^'^-^'^]*  "  "^ol^^^*^'^]-  "  j  2V«     °°  Jo[^o^-(h*h')\] 


(27) 


M-^oFo^^-^'l]*  «^o(^o(^*h')^    - 


-2k  ok 
2ae       ° 


2k  ah 


-2k„ah     Z'    °   .  . 

-  2ae       °  p^  •^o^l^'^^    ^     • 


i'he  last  integral  is  fairly  easy  to  compute   numerically  for  various  values  of 
a  and  p  -  2k  ah.      It  was  found  that  the  values  of  I(h, h)   so  obtained  are  very 
closely  approximated  by  the  tt/2  times  the  value  of  the  integrand  at  n/U. 

In  Figures  1,   2,  3,   U  we  have  graphed  far  various  values  of  a  the  functions 

-p     .» ; 

da 


and 


Here 


F(p)  -  1  *   J  (|)  -  ^2lL    .  2e-P    f  e^J  (i) 
o  a  / — K-  I         o  a 

/l+a  -b 


gO)  .  1    [^Jll^    cos  (±-]  -  i-4    3in(-L)] 


p  -  2k^ah,  F(p)  -  ^  I(h,h) 


0 


and  2+ng(p)  is  the  value  of  the  integrand  in  (26)  for  6  »  n/U.     By  our  approximation, 
WB  should  have 

F(p)     -     1  +  J     g(P). 

The  graphs  show  that  this  relation  is  satisfied  very  closely. 
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When  this  approximation  is  used  in  (25),  the  radiated  pover  per  unit  area  of 
wave front  becomes 


(I6iik  coe)'-'-[2  ■•■  2  cos  2  Jk^  sin  h  ■»•  y>  1 


where 


^^  Y  ■  STK-ITIT  ■  ^  °- 


Note  that  y  is  the  phase  change  on  reflection  from  the  interface.     From  this  formula 
we   see  that  the  radiated  power  will  be  a  minimum  if 


k     sin  h  +  Y 


T,/2 


i.e.,   if 


(28) 


k  h 
o 


/7    cot'"*"  /Z  a. 


For  a  -  1/2,  this  formula  gives  k  h  -  l,35,which  agrees  fairly  well  with  the 
value  1.26     found  by  Cullen.     The  following  table  compares  the  values  for  k  h 
found  by  using  (28)  with  those  found  by  numerical  integration  of   (27). 


Table 

a 

0.2 

O.U 

0.6 

0.3 

1.0 

k  h  from  (27) 

0 

1.62 

i.37 

1.13 

.9U 

0.82 

k  h  from  (28) 

0 

1.33 

1.U7 

1.23 

1.02 

0.88 

We  conclude  that  for  values  of  a  >  O.U,  the  value  h  of  the  most  efficient 
launching  height  as  computed  from  the  approximate  formula  (28)  is  within  10/6  of 
the  height  as  computed  exactly  from  (27), 
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ATTN:     Code  Aer-lD-hlU 

Chief,    Bureau  of  Ships 

Navy  Department 

Washington  25,   D.C. 

AT'm:  ?ir.  E,  Johnston,  Code  838 

Commander 

U.S.  Naval  Air  Missile  Test  Center 

Point  t^gu,  California 

ATTN:      Code  366;   CT-6 


U.S.   Naval  Ordnance   Laboratory 
White  Oak 

Silver  Spring  19,   Maryland 
ATTO:      The  Library 

Commander 

U.S.   Naval   Ordnance 
Test  Station   (Code  753) 
China   Lake,   California 

Librarian 

U.S.  Naval  Postgraduate  School 

Monterey,  Califomia 

Director 

U.S.   Naval  Research  Laboratory 
Washington  25,   D.C. 
ATTO:      Code  1110 

Director,  National  Security  Agency 
Washington  25,   D.C. 

Director 

U.S.   Naval  Research  Laboratory 
Washington  25,  D.C. 
ATTN:     Code  2021 

Director 

U.S.   I^aval  Research  Laboratory 

Washington  25,  D.C. 

ATHJ:     Code  5250,   J.   Bohnert 

Commanding  Officer  and  Director 
U.S.  Navy  Underwater  Sound  Lab. 
Fcrt  Trumbull,  New  London,  Conn- 
ATTN:      Mr.   CM,  Dunn 

Chief  of  Naval  Research 
Navy  Department 
Washington  25,   CC. 
ATTN:      Code  142? 

(2)Commanding  Officer  and  Director 
U.S.  Navy  Electronics  Lab.  lAh, 
San  Diego  52,   California 

Chief,   Bureau  of  Ordnance 
Navy  Department 
Washington  25,   D.C. 
ATTN:      Technical  Library,   Code  Ad3 
Kir.    H.B.  Rex,   Code  ReSh 

Chief,  Bureau  of  Ordnance 

Navy  Department 

Washington  25,  D.C. 

ATTO:   Mr.   C.H.    Jackson,   Code  Re9a 

Chief,   Bureau  of  Ordnance 

Navy  Department 

Washington  25,  D.C. 

ATTN:    'L.E,  Wallingford,    Code  Rehf 

Chief,  Bureau   of  Aeronautics 
Department  of  the   Navy 
Wasliington  25,  D.C. 
ATIW:      Code  Aer-EL-931 

Chief  Bureau  of  Ships 

Navy  Department 

Washington  25,   D.C. 

ATTN:    L.E.   Shoemaker,   Code  838D 

(2) Commander 

U.S.   Naval  Air  Development  Center 
Johnsville,   Pennsylvania 

CoFu-ander 

U.S.   Naval  Air  Test  Center 
Patuxent  River,  Maryland 
ATM:   ET-315,  Antenna   Section 

Naval  Ordnance   Laboratory 
Corona,  California 
ATTN:     Library 

U.S.   Naval  Ordnance  Laboratory 
Department  of  the  Navy 
Corona,   California 
ATTN:     Dr.   Robert  L.  Conger 

U.S.  Naval  Ordnance  Plant 
Indianapolis  18,    Indiana 
ATTI^:     AD-130     (Library) 

Airborne  Instruments  lab..   Inc. 
160  Old  Country  Road 
Mineola,  New  York 
ATTN:     Librarian 
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ACF  Electronics 
800  N.  Pitt  Street 
Alexandria,  Pennsylvania 
ATTN:  Tech.  Library 

American  Machine  and  Foundry  Company 

Electronics  Division 

1085  Commonwealth  Avenue 

Boston  15,  (fess. 

ATTN:     Technical  Library 

Mrs.  Rita  Moravcsik,    Librarian 

Andrew  Alford  Consulting  Engineers 
299  Atlantic  Avenue 
Boston  10,  Massachusetts 
ATTO:     Dr.  A.  Alford 

Battelle  Memorial  Institute 
505  King  Avenue 
Columbus  1,  Ohio 
ATTN:     1-205  F  -K 

Dr.  Bernard  H.  List 

Bell  Aircraft  Corporation 
P.O.  Box  1 
Buffalo  5,  New  York 
ATTN:  Teofmical  Library 

Mrs.  Jasmine  H.  Mulcahey 

Bell  Telephone  Labs.,  Inc. 

Whippany,  New  Jersey 

ATM:     Technical  Information  Lib. 

Bendix  Aviation  Corporation 
Pacific  Division 
11600  Sherman  Way 
North  Hollywood,  Calif. 
ATTN:     Technical  Library 

J.R.   Breninger,  Engr.  Librarian 

Bendix  Radio 

Div.   of  Bendix  Aviation  Corj^oration 
East  Joppa  Road 
Towson  li,   Maryland 

ATTN:     Res.  and  Development  Lab.    (Jj65) 
Chief  Engineer 

Bjorksten  Research  Laboratories,   Inc. 

P.O.  Box  265 

Madison  1,  Wisconsin 

Attni     F.B.  Korsgard,   Treasurer 

Boein?  Airplane  Company 
Pilotless  Aircraft  Division 
Seattle  Ih,  Washington 
ATTN:     Pilotless  Aircraft  Dlv.   Library 
R.R.   Barber,  Library  Supervisor 

Boeing  Airplane  Company 

Wichita  Division 

Wichita,  Kansas 

ATTO:     Kenneth  C.  Knight,   Librarian 

Boeing  Airplane  Company 
Seattle  Division 
Seattle  21i,  Washington 
ATTN:     E.T.  Allen,   Librarian 

Chance   Vought  Aircraft,  Inc. 

P.O.  Box  5907 

Dallas,  Texas 

ATTN:  Engineering  Library 

Mr.  H.S.  White,  Librarian 

Convair,  A  Division  of  Jeneral  Dynamics 

Corp. 
Fort  Worth  1,  Texas 
ATTN:  K.G.  Brown,  Div.  Res.  Librarian 

Convair,  Division  of  General 

Dynamics  Corporat'on 
San  Diego  12,  California 
ATTN:     Engineerin'^  Library 

Cornell  Aero.  Laboratory,   Inc. 

Iih55   Genesee  Street 

Buffalo  21,  New  York 

ATTN:     Elma  T.  Evans,  librarian 

Dalrao  Victor  Company 

IUII4  El  Camino  Real 

San  Carlos,   California 

ATTN:  Dorothea  Qulgley,  Tech.   Libr. 

Dome  and  Mar^lin,  Inc. 
30  Sylvester  Street 
Westbury,  L.I.,  New  York 

Douglas  Aircraft  Company,   Inc. 

Long  Beach  Division 

Long  Beach  1,  California 

ATTN:     Engineering  Library  (C-250) 


Douglas  Aircraft  Comf>any,    Inc. 
B27  Lapha??,  Street 
El  Segundo,  California 
ATTTJ:     Engineerin,:  Library 

Douglas  Aircraft  Company,    Inc. 
3000  Ocean  Park  Bojlevard 
Santa  Monica,   California 
ATTO:     Eq.   Sao.  Reference  Files, 
Eq.  Eng.  A250 

Douglas  Aircraft  Corporation 

Tulsa  Oklahoma 

ATTN:     Technical  Library 

Electronic  Defense  Laboratory 
P.O.   Box  205 

Mountain   View,  California 
ATIN:      Library 

Electronics  Research,   Inc. 
2300  N.  New  York  Avenue 
Evansville,   Indiana 
ATM:     B.H.  Baldridge, 

Engineering  ^pervisor 

Emerson  and  Cuming,   Inc. 
869  Washington  Street 
Canton,  Mass. 
ATTN:     W.R.   Cuming 

The  ^merson  Electric  Mfg.   Company 
8100  Florissant  Avenue 
St.   Louis  21,   I-lissoiiTi 
ATffl:     I-ir.  E.  Breslin, 

Emerson  Engineering  Library 

Falrchild  Aircraft 

Division  Fairchild  Eng.  and  Airplane 

Corp. 

Hagerstown,   Maryland 

ATTO:   Library 


f^amsworth  Electronics  Cr 
Port  Wayne  1,   Indiana 
Arm:     Technical  Library 
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Federal   Telecomnmnication  Labs. 
500  V/ashington  Avenue 
Nutley  10,  New  Jersey 
ATTN:"    Technical  Library 

Gabriel  Electronics 

Division  Gabriel  Company 

Needham  Heights,  Mass. 

ATTN:  Miss  Elaine  Foss,  Librarian 

General  Electric  Company 
Ithaca,  New  York 

ATTN:  Advanced  Electronics  Center 
Library,  John  B.  Travis 

General  Electric  Company 

Electronics  Park 

Syracuse,  New  York 

ATM:  Mrs.  Anne  Damon,  Documents 

Librarian,  Documents  Lib. 

Bldg.  3-11j3 

General  Precision  laboratory.   Inc. 

63  Bedford  Road 

Pleasantville,  New  York 

ATTN:     1-ics,  Mary  0.  Herbst,   Ubrarian 

Goodyear  Aircraft  Corp. 

1210  Massillon  Road 

Akron  15,   Ohio 

ATTN:     Library  DA20  Plant  A 

Grumman  Aircraft  Engineering  Corp. 
Bethpage,   L.I.,  New  York 
ATM:     Engineering  Library, 
Plant  No.  5 

Hallicraf ters  Corporation 
liiiO  West  5th  Avenue 
Chicago,   Illinois 
ATTN:     Library 

Hoffhian  Laboratories,   Inc. 
Los  Angeles,   California 
ATTN :     Librar'.- 


Hughes  Aircraft  Company 
Microwave  Laboratory 
Culver  City,  California 
ATTO:    Library,   Building  12 

M.D.  Adcock,   Head,   Antenna  Dept. 

Hughes  Aircraft  Company 
Culver  City,  California 
ATTN:  Research  and  Development  Lib. 

ftycon  Eastern,   Inc. 

75  Cambridge  Parkway 

Cambridge,  I-Iass. 

ATTN:     Mr.   3.   Stein,   Library 

International  Business  Hacines  Corp. 
I-Iilitary  Products  Division 
590  i'iadison  Avenue 
New  York  22,  liew  Y-^rk 
ATTN:     Mr.   C.F.  McElwain 

International  Business  Machines  Corp. 
Military  Products  Division 
Airborne  Computer  Laboratories 
/estal.   New  York 
ATM:     Mr.   J.J.  Walsh,  Dept.    528 

Jansky  and  Bailey,    Inc. 

1339  Wisconsin  Avenue,  N.W. 

Washington,  D.C. 

ATTN:   Tech.   Library,  Delmer  C.   Ports 

Henry  Jasik,  Consultin,'.  Engineering 
61  Second  Street 
I'dneola,  Sew  York 

Jet  Propulsion  Laboratorv 
U800  Oak  Grove  Drive 
Pasadena,  California 
ATTN:      I.E.   Newlan 

Martin  Katzin  and  Co.,   Consultants 
711     ll4th  Street,  N.W. 
Washington  5,  D.C. 

Lincoln  Laboratory 

P.O.   Box  73 

Lexington  73,   Massaohuset's 

ATTO:   Henry  Straus,  Document  Rm.,  A-229 

Lockheed  Aircraft  Corporation 
f-bssiie  Systems  Division 
7701  Woodley  Avenue 
Van  Nuys,  California 
ATTN:  MSD  Library  -902,   Dept.   82-10 
E.A.  Blasi,  Head,  Antenna  Ub. 

!<3c.^heed  Aircraft  Corporation 

California  Division 

Engineering  Library 

Dept.  72-25,   Plant  A-1,   Bldg.   63-1 

ATTN:     N.C.   Harnois 

The  Glenn  L.  Martin  Company 

Denver,  Colorado 

ATTN;  Wm.   V.    Foley,  Engineering  Lib. 

The  Glenn  L.   Martin  Company 
Baltimore  3,   t-laryland 
ATTO:  Engineering  Library 
Antenna  Design  Group 

Maryland  Electronic  Mfg.  Corp. 
5009  Calvert  Road 
College  Park,   I-'aryland 
ATIT):     H.  Warren  Cooper 

Mathematical  Reviews 
80  Waterman  Street 
Providence  6,  Rhode  Island 

The  W.L.  Maxon  Corporation 
Ii60  West  31ith  Street 
New  York,  New  York 
ATTN:     Antenna  Laboratory 

McDonnell  Aircraft  Corporation 

Box  516 

St.  Louis  3,  Missouri 

ATTN:     R.D.  Cetrich 

Engineering  Library 

teMillan  Laboratory,   Inc. 
Ipswich,   Massachusetts 
ATTN:     Security  Officer 
Docume'^t  Room 

Melpar,  Inc. 

3000  Arlington  Boulevard 

Falls  Church,   Virginia 

ATTN:  Engineering,  Technical  Lib. 


Micrbwave  Radiation  Company 
19223  So.    Hamilton  Street 
Gardena,   California 
ATTN:      Mr.    Morris   J,  Ehrlich,    Fres. 

Northrop  Aircraft,   Inc. 
Engineerinj^  Department 
lOCC  E.  Broadway 
Hawthorne,   California 
ATTN:     E.A,  Freitas,  Engineering^ 
Library,   Plant  III 

North  American  Aviation,   Inc. 
12?ll4  Lakewood  Boulevard 
Downey,   California 
ATTIh      Technical   Library, 
Dept.  991-h 

North  American  Aviation,    Inc. 
Los  Angeles   International  Airport 
Los  Angeles  hS,   California 
ATTN:     Engineering  Technical  File 

Page  Communications  Engineers 
71c  Ihth  Street,   N.W. 

VJashlnpton   5,  D,G, 
ATTN:      Library 

Phllco  Corporation  Research  Division 
A  and  Lippincott  Street 
Philadelphia  3lj,  Pa. 
ATTO:     Technical  Library 

Pickard  and  Bums,   Inc. 

2ljO  Highland  Avenue 

Needham  9h,  Mass. 

ATTHz     Dr.   J.T.   de  Bettencotu-t 

Polytechnic  Research  and  Development 

Con!pany,    Inc. 

2C2  Tiilary  Street 

Brooklyn  1,  New  York 

ATIW:     Technical  Library 

Radiation  En^neering  Laboratory 
Kaynard,   I-Iass. 
ATTN:     Dr.   John  Ruze 

Radiation,    Inc. 
P.O.  Drawer  37 
Melbourne,  Florida 
ATTN:      Technical  Library 

Mr.  M.L.   Cox,  Assistant 

Project  Engineer 

Piadio  Corp.   of  America 

RCA  Laboratories 

Rocky  Point,   New  York 

ATTN:     P.S,  Carter,  Lab.   Library 

RCA  Laboratories 
Princeton,   New   Jersey 
ATTN:      Miss   Fern  Cloak,   Librarian 
Research  Library 

RCA  Defense  Electronic   Products 
Canden  2,   New  Jersey 
ATTN:     Miss   J.M.   Steever,   Librarian 
Library  10-2 

(2jrhe  RajTO-^ooldrid^e  Corporation 
S73C  Arbor  Vitae  Street 
Los  Angeles  U5,   California 
ATTN:      I-^r^aret  C,  Whitnah, 

Chief  Librarian,  Research 
and  De\'elopment  Library 

The  Rand  Corporation 

Via:      3BAMA  Liaison  Office 

1700  Main  Strer.t 

Santa  Monica,  California 

ATTN:      Library 

Rantic  Corporation 
Carlabasos,  California 
ATTN:  Technical  Library 

Raytheon  Manufacturirij^  Company 
M  and  R  Division 
Bedford,  Massachusetts 
ATTW:      Irving  Goldstein 

Rayth.eon  Laboratory 
Wayland,  Massachusetts 
ATIT'I:     Document  Lib.,    John  E.  Walsh 
Dept.   897U,  Antenna 
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Raytheon  Manufa -turing  Company 
Wayland  Laboratory 
Wayland,   Massachusetts 
ATTN:      Miss  Alice  G.   Anderson,   Lib. 
Wayland  Library 

Republic  Aviation  Corporation 

Farmingdale,  L.I.,  New  York 

ATTTJ:   H.  Stleglltz,  Engineering  Lib. 

Rheem  Manufacturi  ne'  Company 
Research  and  Development 
9236  East  Hall  Road 
Downey,  California 
ATTN:   J.C.  Joerger 

Ryan  Aeronautical  Company 
San  Diego,  California 
ATTN:   Technical  Library 

Sanders  Associates 
137  Canal  Street 
Nashua,   New  Pianipshire 
ATTN:     N.R.  Wild,   Library 

Sandia  Corporation 

P.O.  Box  5800 

Classified  Document  Division 

Albuquerque,   New  Mexico 

Sperry  Gyroscope  Company 
Great  Neck,   L.I.,  New  Yoiic 
ATTN:     Engineering  Library 

Florence  W,   Tumbull, 

Engr.   Librarian 

Stanford  Research  Institute 
Menlo  Park,   California 
ATTTJ:     Engineering  Document  Library 
Mrs.   Beatrice  Gibson,   Libr. 

Sylvania  Electric  Company 
100  First  Street 
Waltham,   Mass. 
ATT:      Librarian 

Systems   Laboratories  Corporation 
I50I6    Ventura  Blvd. 
Sherman  Oaks,  California 
ATITJ:     Dr.  W.C.  Hoffman 

Technical  Research  uroup 
17  Union  Square  West 
New  York  3,  New  York 

A.S,   Thomas,    Inc. 
161  Devonshire  Street 
Boston,  Mass. 
ATTN:     A.S.   Thomas 

Westlnghouse  Electric  Corp. 
2519  Wllkens  Avenue 
Baltimore  3|   Maryland 
ATTN:     Engineering  Library 
Kent  M.   Mack 

Wheeler  laboratories.   Inc. 
122  Cutter  mil  Road 
Great  Neck,   New  York 
ATTN:     Mr.   Harold  A,  Wheeler 

Zenith  Plastics  Comoany 

rOX    91 

Gardena,  California 
ATTN:     Mr.   S.3.   Oleesky 

Library 

Geophysical   Institute 
University  of  Alaska 
College,  Alaska 

Electronic  Research  Laboratory 
Div.   of  Electronic  Engineering 
University  of  California 
Berkeley  h,   California 
ATra:     S.  Silver 

Electronics  Research  Laboratory 
232  Cory  Hall 
University  of  California 
Berkeley  h,   California 
ATTI^:      J.R.  Whinnery 

California  Institute  of  Technology 
1201  E,  California  Street 
Pasadena,  California 
ATTiit     Dr.   Charles   H.   Papas 


Carnegie  Institute  of  Technology 
Pittsburgh  13,   Pennsylvania 
ATTN:     Albert  E,   Heins 

School  of  Electrical  Engineering 

Cornell  University 

Ithaca,   New  York 

ATHJ:     Dr.   Henry  G.   Booker 

College  of  Engineering 
University  of  i-l.orida 
Gainesville,   Florida 
ATTN:     Engineering  Sciences  Library 
Professor  M.H.   Latour 

Engineering  Experiment  Station 
Georgia  Institute  of  Technology 
Atlanta,  Georgia 
ATITl!       Georgia  Tech.  Library 
Mrs.    J.H.   Crosland 

Technical  Reports   Collection 

Har-.-ard  University 

Gordon  McKay  Library 

303A  Pierce  Hall 

Oxford  Street 

Cambridge  38,   Massachusetts 

ATTN:     Mrs.  E.L.   Hufschmidt,   Librarian 

Harvard  College  Observatory 

60  Garden  Street 
Cambridge,   Massac hiisetts 
ATTN:     Dr.   Fred  L.  Whipple 

Serials  Department     -220-S  Library 

University  of  Illinois 

Urbana,   Illinois 

ATTN:      Engineering  Library 

Antenna  Section 

Electrical  Engr.  Research  laboratory 
University  of  Illinois 
Urbana,   Illinois 
ATTW:     EERL  Library,   Room  21C, 
Dr.   R.H.  DuHamel 

The  Johns  fiopkins  University 
Department  of  Physics 
Baltimore  16,  Maryland 
ATTN:     Professor  Donald  E.  Kerr 

Radiation  Laboratory 
Johns  Hopkins  University 
131s  St.   Paul  Street 
Baltimore  2,  I-laryland 

Applied  Physics  laboratory 
The  Johns   Hopkins  University 
8621  Georgia  Avenue 
Silver  Spring,  I-laryland 
ATTN:     Document  Library 

Massachusetts  Institute  of  Technology 
Research  Laboratory  of  Electronics 
Document  Room  20B-221 
Cajnbridi?e  39,   Massachusetts 

University  of  itLchigan 
Electronic  iJefense   li-oup 
Engineering^  Research  Institute 
Ann  Arbor,  Michigan 
ATTN:      J. A,   Boyd,   Supervisor 

Willow  Run  Labs.,   of  the  Engineering 
Research  Institute 
The  University  of  I-lichigan 
Willow  Run  Airport 
Ypsllanti,  Miichigan 
ATTO:      K.M.   Sie;?el,   Head,   Theory 
and  Analysis  Department 

University  of  Michigan 
Willow  Run  Laboratories 
Willow  R\in  Airport 
Tpsilanti,   Michigan 
ATTN:     Sheila  Coon,   Librarian 

University  of  Minnesota 
Minneapolis ,   Minnesota 
Engineering  Library 
ATTN:      HeaH  Lihrarian 
Microwave  Laboratory 
Electrical  Engineering  Dept. 
Northwestern  University 
Evanston,    Illinois 
ATTN:     Prof.  H.E,  Beam 
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Antenna  Laboratory 

Department  of  iilectrica]   Engineering 

Ohio  State  University 

Columbus,  Ohio 

ATTN:     Dr.   T.E.  Tice 

Oklahoma  University  Research  Found. 

Norman,  Oklahoma 

ATTN)  Technical  Library 

Professor  C.L.  Farrar 

Microwave  Research  Institute 
$5  Johnson  Street 
Brooklyn  1,  New  Y^rk 
ATOJ:  "Dr>  A. A.   Oliner 

Polytechnic  Institute  of  orooklyn 

55  Johnson  Street 

Brooklyn  1,  New  York 

ATITJ:   Microwave  Research  Institute 

Syracuse  University 
Dept.  of  Electrical  Engr. 
Syracuse  10,  New  York 
ATTN:     Technical  Library 

University  of  Texas 
Electrical  Engr.  Research  Lab# 
P.O.  Box  8026,  Univ.   Station 
Austin  12,    Ttexas 
ATIN:      John  R.   Gerhardt 

Dept.   of  Physics 
University  of  Texas 
Austin  12,  Texas 
ATTN:     Physics  Library 
Claude  W.   Horton 

Tufts  University 

Medford  55,  Massachusetts 

ATTN:     Research  Ub.   of  Physical 

Electronics 

Professor  Charles  R.  tiinr^ins 

University  of  Toronto 
Dept.   of  Electrical  Engineering 
Toronto,   Ontario,   Canada 
ATTO:   Prof.   'George  Sinclair 

University  of  Washington 
Dept.   of  Electrical  Engineering 
Seattle,  Washington 
ATTN:      Prof.    3.    Held 

Ionosphere  Research  Laboratory 
Pennsylvania  State  College 
State  College,   Pennsylvania 

ATTN:   Prof.  A.H.  Waynick,  Director 

Institute  of  MatheTTatir.i.1   Science 

25  Waverly  Place 

New  York  3,  New  York 

ATTN:     I-Irs.   JoAn  Segal,   Librarian 

School  of  Electrical  Engineering 

Purdue  University 

lafayette,  Indiana 

ATTN:     Professor  F.V.  Schultz 

Department  of  Electrical  Engineering 
California  InsUtute  of  Technology 
Ai'TN:  Dr.  Charles  H.   Papas 

Electronics  Division 
Rand  Corj-^oration 
1700  Main  Street 
Santa  Monica,  California 
ATTO:     Dr.  Robert  Kalaba 

National  Bureau  of  St*  dards 

Washington,  D.C. 

ATTN:  Dr.  W.K.  Saunders 


Applied  Mathematics  and  Statistics 

Stanford  University 

Stanford,  California 

AWN:  Dr.  Albert  H.  Bowker 

Department  of  Physics  and  Astronomy 
Michigan  State  University 
East  Lansing,  Michigan 
ATTO:     Dr.  A.  Leitner 
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University  of  Tennessee 
Knoxville,  Tennessee 
ATTN:  Dr.  Fred  A,  flcken 

California  Institute  of  Technology 
1201  E.  California  Street 
Pasadena,  California 
ATTN:  Dr.  A.  Erdelyl 

Wayne  University 

Detroit,  Michigan 

ATTN:  Professor  A.F.  Stevenson 

Mathematics  Department 
Stanford  University 
Stanford,  California 
ATTN:  Dr.  Harold  Levine 

University  of  Minnesota 
Minneapolis  lli,  Minnesota 
ATM:  Prof.  Paul  C.  Rosenbloom 

Department  of  MatJiematics 
University  of  Pennsylvania 
Philadelphia  h.   Pa. 
ATTN:  Professor  Bernard  Epstein 

Applied  Physics  laboratory 
The  Johns  Hopkins  University 
8621  Georgia  Avenue 
Silver  Spring,  Maryland 
ATTN:  Dr.  B.S.  Oouray 

(2>Exchange  and  Gift  Division 
The  Library  of  Congress 
Washington  25,  D.C. 

Electrical  Engineering  Department 
Massachusetts  Institute  of  Tech. 
Cambridge  39,  f-lassachusetts 
ATTN:  Dr.  L.J.  Chu 

Nuclear  Development  Assoc.  Inc. 

5  New  Street 

White  Plains,   New  York 

ATTN:      Library 

Lebanon  Valley  College 
Annville,  Pennsylvania 
ATTN:     Professor  B.H.    BissiniTer 

Dept.  of  Physics,   Thaw  Hall 
University  of  Pittsburgh 
Pittsburgh  13,   Pa. 
ATTN:     Dr.  Edward  Gerjuoy 

Dept.   of  Physics 
Amherst  College 
Amherst,  Massachusetts 
ATTN:     Dr.  Arnold  Arons 

California  Institute  of  Tech. 
Electrical  Engineering 
Pasadena,   California 
ATM:     Dr.   Zohrab  A.  Kaprlellan 
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